Role of Particle Interactions in the Feshbach Conversion of Fermion Atoms to Bosonic 

Molecules 
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We investigate the Feshbach conversion of fermion atomic pairs to condensed boson molecules 
with a microscopic model that accounts the repulsive interactions among all the particles involved. 
We find that the conversion efficiency is enhanced by the interaction between boson molecules while 
suppressed by the interactions between fermion atoms and between atom and molecule. In certain 
cases, the combined effect of these interactions leads to a ceiling of less than 100% on the conversion 
efficiency even in the adiabatic limit. Our model predicts a non-monotonic dependence of the 
efficiency on mean atomic density. Our theory agrees well with recent experiments on 6 Li and 40 K. 
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the research activities in cold atom physics [J, [2, S 0] 
since its first experimental realization^. Among these 
research activities, the production of diatomic molecules 
from Fermi atoms with Feshbach resonance is of special 
interest and has attracted great attention. First, it is 
an interesting phenomenon by itself; second, it provides 
a unique_ experimental access to the BCS-BEC crossover 
physicsQ. So far, by slowly sweeping the magnetic field 
through the Feshbach resonance, samples of over 10 5 
weakly bound molecules (binding energy ~ 10 kHz) at 
temperatures of a few tens of nK have been produced 
from quantum degenerate Fermi gasp], la, la]- 

The Feshbach conversion is a complicated process in- 
volving many fermion atoms and boson molecules in a 
sweeping magnetic field that crosses a resonance. The 
theoretical description of the conversion efficiency as 
a function of sweep rate, atom mass, atomic density, 
and temperature is still under development. The ex- 
isting theories include the Landau-Zener (LZ) model of 
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two-body molecular productionflC 
body extension at zero temperature 
space density modelfl5j. and equilibration model at finite 
temperatures [l6j . 

In this Letter we study a microscopic model of the 
Feshbach conversion that accounts all the two-body in- 
teractions, which include atom-atom, molecule-molecule, 
and atom-molecule interactions. These interactions are 
ignored in previous theoretical studies [H GJ, EH 0- 
We find that these interactions affect strongly the Fesh- 
bach conversion efficiency: the repulsive interaction be- 
tween molecules tends to enhance the conversion effi- 
ciency while the other two repulsive interactions between 
atoms and between atom and molecule suppress the effi- 
ciency. Combined together, these interactions can yield a 
ceiling of less than 100% for the conversion efficiency even 
in the adiabatic limit of the sweeping magnetic field. This 
interaction-suppressed conversion efficiency is in spirit 
the same as the broken adiabaticity by interaction in the 
nonlinear LZ tunneling [I?], [l8j]- I n addition, our model 



predicts a non-monotonic dependence of the conversion 
efficiency on mean atomic density. Our results are com- 
pared to recent experiments with 6 Li and 40 K[8|.Ii^: they 
are in good agreement. 

To include all particle interactions, we extend the two 
channel model[l|| H(| HH and write the Hamiltonian as 
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Here £k = h 2 k 2 /2m a is the kinetic energy of the atom 
and a —\, J, denote the two hyperfine states of the atom. 
U a = AU , Uab = AU x ,g = Ago, and 7 = 70 - Ag%/U c , 
where 70 = ^ co {B — -Bo) is the molecule energy under 
the linearly changing magnetic field with B = —a r t, 



go = y 4:7rh ab g ABfj, co /m a is the atom-molecule cou- 
pling 3], Uq — 4:Trh 2 cibg I m a is the interaction between 
atoms, U\ = 4:7rh 2 l. 2ab g /m a b is the atom-molecule scat- 
tering interaction, and Ub = 4Trh 2 0.6ab g /mb is the in- 
teraction between molecules 22] . A = (1 + Uq/Uc)^ 1 
and U^ 1 = - Sfc( ex P( — k 2 /k1)/2e k ) with the cutoff mo- 
mentum k c representing the inverse range of interaction 
21, 23, 24]. Bo and AB arc Feshbach resonance point 
and width, respectively. m a and m?, = 2m a are masses 
for atoms and molecules, and m a b = ^m a is the reduced 
mass for the atom-molecule interaction. 

In experiments, the molecular bosons are more tightly 
confined in space than the fermion atoms due to their 
different statistics. To reflect this, we use V a for the 
volume of fermion atoms and Vb for bosonic molecules. 
We assume the zero temperature limit, consider only one 
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bosonic mode, and ignore all possible dissipations in the 
system, such as the loss of atoms by three-body collision. 

In current experiments, the intrinsic energy width 
of a Feshbach resonance is much larger than the 
Fermi energy -Ef[25|, it is therefore reasonable to 
assume ek = e. This is called degenerate model 
in Ref.pJ, 0J, H3]. We introduce the following 

Sk( a k,T a -k,| b + bta -k,4 a k,T) 
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With the commutators [X z , L^. 

can obtain the Heisenberg equations for the system 
ih-j^Li — [Li, H],(l = x, y, z). Since all the commutators 
vanish in the limit of N — > 00 and N is large in current 
experiments, it is appropriate to take L Xl L y , and 



L z as three real numbers u, v, w, respectively. 
Heisenberg equations are then reduced to 



These 



du/dr 



-Sv - 2xv(l - w) , dw/dr = V2v 1 (2) 



3\/2 1 
dv/dr = -t-(w-l)(w + -)+5u + 2xu(l-w),(3) 

where r = (gVb\fN /KVa^ 2 )t. Because of the identity 
u 2 +v 2 = l) 2 (w+l), with introducing the canonical 

variable 9 = arctan(t> /u) we have a classical Hamiltonian, 



H = Sw - x(l - w) 2 + v 7 ^ - !) 2 ( w + 1) co s 0. (4) 

The above equations show that all the experimental pa- 
rameters affect the system via only two dimensionless 
parameters <5 and X - By a trivial shift of time origin, we 
can set 5 = ar with 



a r AnhnabgAB ^ 2 V b 2 
a m n V 2 



(5) 



where n — N/V a is the mean atomic density. The non- 
linear parameter x is given by 



X = 2.3- 




0.l5V a 



(0) 



To understand the dynamics, we first look at the fixed 
points w = u = v = of Eqs. (|2l3p . The energies for these 
fixed points make up energy levels of the system as shown 
in FigQ] One sees that the structure of these energy 
levels changes dramatically as the nonlinear parameter 
X increases. Specifically, we observe: (?) There are two 
fixed points when \S\ is large enough: one for bosonic 
molecule (BM) and the other for fermion atom (FA), (ii) 
When \5\ < 5 C = V2, there is an additional fixed point 
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FIG. 1: Adiabatic energy levels for different interaction 
strengths. (a) X = 0; (b) X = Xc = V2/4; (c) X = 1.5. The 
unstable states are indicated by dashed lines (MQ and DM). 



with w = 1. However, this fixed point is dynamically 
unstable [Hj]. (Hi) For X > Xc — there appears 

one more fixed point denoted by P3 and, consequently, 
a loop in the energy levels. As we shall see, this loop 
has highly non-trivial physical consequences. This fixed 
point P3 is also unstable. 

Consider the adiabatic evolution of the system start- 
ing from a high negative value of S with w = 1. This 
corresponds to the experiments where the magnetic field 
sweeps slowly across the Feshbach resonance with no 
bosonic molecules initially. When X 1S small, such as 
in FigfjTa), the evolution of the system follows the solid 
line, converting all fermion atoms into molecules. How- 
ever, when x is beyond Xc as in Fig[T](c), the system will 
find no stable energy level to follow at singular point M. 
As a result, only a fraction of fermion atoms are con- 
verted into bosonic molecules. 

This simple analysis is confirmed by our numerical 
results, which are plotted in FigO In our calcula- 
tion, the 4-5th Runge-Kutta step-adaptive algorithm is 
used in solving the differential equations (I2|3|) . Because 
w = 1 is a fixed point when 6 < —y/2, we start from 
(id, u, v) ~ (1, 0, 0) and sweep the field from 8 = —\pl to 
200. In Fig[21 the conversion efficiency T, i.e., the frac- 
tion of the converted fermion atom pairs is drawn as a 
function of a. Evidently, T approaches one as a — > 
when x < Xc, indicating that all atomic pairs are con- 
verted into molecules. In contrast, when X > Xc T does 
not increase to one in the adiabatic limit a — * 0. This 
means that there is a ceiling T a d(< 100%) on the conver- 
sion efficiency. Moreover, Fig|2] demonstrates that pos- 
itive x suppresses the conversion efficiency whereas the 
negative x enhances it. Because the repulsive interaction 
between bosonic molecules enters x as a negative value, it 
enhances the conversion efficiency; the repulsive fermion 
atomic interaction and atom-molecule interaction con- 
tribute positively to x, they suppress the conversion. 

The ceiling T ac i on the atom-molecule conversion ef- 
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(a) x=OA S=-U> (b) x =0.6, 5=-^2 (c) ^=0.6, 6=0 
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FIG. 3: Phase spaces of Hamiltonian @. The dark line in (a) 
is for the fixed point w — 1, it = 0, v = 0. It is a line because 
9 is not defined at u — v = 0. The two fixed points on line 
w — 1 in (c) are in fact the same fixed point; they are artifact 
caused by the definition 6 = arctan(u/u). 



FIG. 2: Conversion efficiency F as a function of the sweeping 
rate a for various interactions. 



ficiency depends on \. This dependence can be found 
by examining the phase space diagrams of our system 
shown in Fig|3] As 8 ramps up slowly from a large neg- 
ative value, the fixed point P3 will move up until it hits 
the fixed point w = l,u = 0,v = 0, represented by a 
dark straight line in Fig[3ja). This collision occurs at 
S = —\/2. Immediately after the collision, the hyperbolic 
fixed point P3 is no longer a fixed point and becomes 
a solution that evolves along the dark line in Fig(3£b). 
The dark line is given by v2 = %(1 — w) — y/1 + w cos#, 
which is found by taking E = 5 = —\[2 in the Hamil- 
tonian (HJ). As the action of this trajectory is nonzero 
while a fixed point has zero action, this collision of the 
two fixed points represents a sudden jump in action. It 
is this sudden jump that has caused the nonzero fraction 
of remnant atoms. As 6 ramps up further slowly, the 
trajectory will change its shape as witnessed in FigJ3]Jc); 
however, its action stays constant as demanded by the 
classical adiabatic theorem! 271 12811. The action is 



1 

2tt 



2tt JC 



cos 6 V 8x 2 -4V2x+cos 2 8 
2tt 4 X 2 -2V^ X +cos 2 j q 

2^2 U17 l 



which yields the ceiling on the efficiency 
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Now we compare our theory with existing experi- 
ments. For the experiment with 6 LiQ, the mean den- 
sity is n = 4 x 10 12 cm~ 3 with N — 6 x 10 5 atoms. 
The scattering length cib g = 59a b, (J-co ~ 2/^s, where 
as and \xb are Bohr radius and Bohr magneton, re- 
spectively, and the resonance width AB = 0.1G at 
Bo = 543. 8G. The Fermi energy Ep in the combined 
harmonic and box-like trapping potential of Ref. [8J is 

given by E F = [WirNh 3 ^/ (8y/2m a L)] 2 ^ 5 , where to r = 
2ir x 800s" 1 is the angular frequency of the radial har- 
monic trap and L — 480/im is the size of the axial po- 



tential. The ground state energy of molecular bosons 
is E G = hcu r + h 2 Tr 2 /(2m b L 2 ). Then we have V a /V b = 
Ep/Ec = 36. We set A = 391 with a momentum cut- 
off K c — 96kp 29]. From Eq.©, the sweeping rate is 
a r /a = 20G/ms. The second term in the bracket of 
Eq. ^ that accounts for the repulsive interaction between 
bosonic molecules is small, so the interaction parameter 



takes the form of x 



2 3^ 



From the above 



experimental parameters, we find the interaction param- 
eter as x = 1-26. This strong interaction (> X c) indicates 
a ceiling of T a( i = 0.48 via Eq.©. This is in good agree- 
ment with experiments (see Fig0^). 

For 40 K, the situation is different. The resonance at 
B = 202. 1G has a large width of AB = 7.8G and 
the mass of 40 K is 7 times that of 6 Li. In Ref. (l5| . 
the fermions are confined in a dipole trap character- 
ized by a radial frequencies v r between 312 and 630 
Hz and an aspect ratio of v r jv z — 70. The Fermi en- 

1/3 

ergy is Ep — Ti (3Nuj 2 lj z ) and the ground state en- 
ergy of condensed bosons is Eq = huj r + hu z /2. For 

the dipole trap, the ratio V a /V b = {E F /E G ) 3/2 = 251. 
With a bg = 174as, /i co ~ 2/j,b, initial clouds have mean 
densities n = 2 x 10 12 cm -3 , and N = 2.5 x 10 5 [l~5l|, 
we obtain X — 0.135 ~ 0.274, which is less than the 
threshold Xc = V2/i. Therefore, 40 K atom pairs can be 
completely converted to bosonic molecules in adiabatic 
limit. Indeed, the conversion efficiency up to 90% has 
been observed [l5l]. 

We emphasize that the suppressed conversion efficiency 
by particle interaction dominates only at low tempera- 
tures. As a result, in the above we have only compared 
to the data obtained at low temperatures (T/Tp = 0.1 
for 6 Li and T/T F = 0.05 for 40 AT). Temperature can af- 
fect the conversion efficiency strongly as reported in Ref. 
[lo] ] . The ceiling of 50% conversion efficiency observed in 
Ref. [tJ is likely a thermal effect since the experiment is 
performed at TjTp — 0.33, and has been explained by 
the theories of finite temperature [30l. [3l| 

For the 6 Li, from Eq.(2)(3)we have also calculated nu- 
merically the conversion efficiency as a function of sweep- 
ing rate. The comparison between our theory and exper- 
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FIG. 4: (a) Comparison between our theory and experimental 
data of 6 Li;S] for the conversion efficiency T as a function 
of field sweep rates, (b) The dependence of V on the mean 
atomic density. 



iment is shown in Fig[4^,. They are in a good agreement. 
In addition, our model predicts a non-monotonic depen- 



dence of the conversion rate on the mean atomic density 
(see FigBJb) ) . This can be understood from Eqs . (|5l6p . In 
Eq.©, we see the effective sweeping rate a is inversely 
proportional to the atomic density. So, increasing the 
density will reduce the effective sweeping rate and there- 
fore enhance the conversion rate. On the other hand, 
higher density will give larger nonlincarity \ as indicated 
in Eq.([5]), which in turn suppresses the atom- molecule 
conversion. These two factors compete with each other, 
giving rise to the non-monotonic curves in Fig[4fb) . In 
practical experiments, to achieve higher conversion effi- 
ciency, one needs to carefully choose initial fermion atom 
density, making it fall into the optimal parameter regime. 

In summary, we have identified the significant role of 
the interactions between particles in the Feshbach con- 
version of atomic fermion pairs to molecular bosons. Our 
theory is consistent with the existing experiments. Our 
model also predicts a non-monotonic dependence of the 
conversion rate on the mean atomic density, which is im- 
portant for the optimal choice of parameters in future 
Feshbach experiments. 
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